New cytotoxic dammarane type saponins from Ziziphus spina-christi

Cancer is the world's second-leading cause of death. Drug development efforts frequently focus on medicinal plants since they are a valuable source of anticancer medications. A phytochemical investigation of the edible Ziziphus spina-christi (F. Rhamnaceae) leaf extract afforded two new dammarane type saponins identified as christinin E and F (1, 2), along with the known compound christinin A (3). Different cancer cell lines, such as lung cancer (A549), glioblastoma (U87), breast cancer (MDA-MB-231), and colorectal carcinoma (CT-26) cell lines, were used to investigate the extracted compounds' cytotoxic properties. Our findings showed significant effects on all the tested cell lines at varying concentrations (1, 5, 10, and 20 µg/mL). The three compounds exhibited potent activity at low concentrations (< 10 μg/mL), as evidenced by their low IC50 values. To further investigate the complex relationships between these identified cancer-relevant biological targets and to identify critical targets in the pathogenesis of the disease, we turned to network pharmacology and in silico-based investigations. Following this, in silico-based analysis (e.g., inverse docking, ΔG calculation, and molecular dynamics simulation) was performed on the structures of the isolated compounds to identify additional potential targets for these compounds and their likely interactions with various signalling pathways relevant to this disease. Based on our findings, Z. spina-christi's compounds showed promise as potential anti-cancer therapeutic leads in the future.


Extraction and fractionation of Ziziphus spina-christi
Fresh Z. spina-christi leaves (1 kg) were harvested and air-dried -for one week in the shade.The finely powdered leaves were macerated in 70% methanol (5 L, 3 days each) at room temperature, and the crude extract was concentrated using a rotary -evaporator (Buchi Rotavapor -R-300, Cole-Parmer, -Vernon Hills, -IL, USA) under vacuum at 45 °C.The dried extract was suspended in 500 mL of distilled water (H 2 O), then portioned with various polarity solvents (n-Hex., -DCM, EtOAC, -and n-but.).In each step, the organic phase was evaporated at decreased pressure to obtain fractions I (10.0 g), II (5.0 g), III (12.0 g), and IV (20.0 g), while the -remaining mother -liquor was concentrated to yield the -aqueous fraction (V).For biological and phytochemical research, all fractions were kept at 4 °C46-52 .

Cell viability testing using MTS
A549, U87, MDA-MB-231, and CT-26 cancer cells were plated in 96-well plates at a density of 2X10 4 cells/well in 100 μL of medium before the addition of different concentrations of the plant extracts.Untreated control group was incubated with -100 µL/well of fresh media and -methanol.After 24 h., the medium was aspirated and replaced by -100 µL of fresh media, and 20 µL of MTS -reagent in each well (CellTiter 96 -Aqueous One -Solution cell -proliferation assay, Promega -Corporation, Madison, -WI, -USA).The plates were then incubated at 37 °C with 5% CO 2 for 2 h.The absorbance was measured at 490 nm using a Spectramax plus 384 Microplate reader (Molecular Devices, Sunnyvale, CA, USA).Relative cell viability values were expressed as the percentage of the absorbance from the treated wells compared to the control wells (untreated), given that the control wells viability was set to 100%.The cells were examined under a cell imaging system (EVOS FL Digital Microscope) using a 20X objective as well before getting rid of the plates.The half-maximal inhibitory concentrations (IC 50 ) of pure compounds for four human cancer cell lines were calculated.

Molecular docking MD simulation and network construction
Docking was performed using Aut -Dock Vina -software, MD -simulations were performed -using NAMD 3.0.0software, and PPI network building was performed using cytoscape.The companion file contains detailed descriptions of these processes.

Phytochemical investigation of Ziziphus spina-christi
Sugar moieties in saponins 1, 2, and 3, were identified by GC-MS after acid hydrolysis followed by derivatization with L-cystein methyl ester, subsequent silylation, and comparison with derivatized reference sugars.
The positive ion HRESIMS spectrum of

Cytotoxic activity
The ability to make soap distinguishes steroidal or triterpenoid glycosides as saponins.Different saponins have been identified and purified, and their use in cancer chemotherapy is growing.High structural diversity in saponins is associated with anticancer effects such as reversing multidrug resistance (MDR), anti-metastasis, anti-proliferation, and anti-angiogenesis.These effects are brought about by the promotion of cell differentiation, bile acid-binding, induction of apoptosis, immune-modulatory effects, and amelioration of carcinogen-induced cell proliferation 60 .Consequently, different molecular mechanisms are contribute to saponins' anticancer action.It should be noted that the aglycone moiety, the length and linkage of the glycosidic chain, the presence of a functional carboxylic group on the aglycone chain, the number of sugar molecules and hydroxyl group, the position of the hydroxyl group, stereo-selectivity, and the type of sugar molecule on the glycine chain are all structural components that have a significant impact on the mechanism of saponins' anticancer action [61][62][63] .
In this investigation, the isolated dammarane type saponins compounds 1-3 from the methanolic extract of Ziziphus spina-christi leaves were evaluated for their cytotoxic activities against A549, U87, MDA-MB-231, and www.nature.com/scientificreports/CT-26 cancer cell lines, using MTS.The results showed that three pure compounds exhibit a high anticancer efficacy, as evidenced by their ability to significantly reduce cancer cell viability.Notably, these compounds were strongly active at concentrations below 10 μg/mL, resulting in a 50% cytotoxic effect (Table 3, Fig. 3).Dammarane-type saponins gained attention after multiple articles demonstrated their anticancer activity against a range of malignancies.Gypenosides CP2, 5, 6, 7, 11, 12, 13, and 15 were identified and characterised from Gynostemma pentaphyllum aerial parts and demonstrated antiproliferative activities against two human tumour cell lines (A549 and HepG2) 64 .The EC 50 values of these eight compounds against HepG2 cells ranged from 29.3 to 100.6 μM, while only four compounds (2, 11, 13, and 15) had EC 50 values less than 100 μM (EC 50 59.487.3 μM) against A549 cells.Compound 2 was the most effective against A549 cells, whereas Compound 5 was among the most effective against HepG2 cells 64 .Another dammarane-type saponin, gypenoside Jh1, isolated from heat-processed G. pentaphyllum ethanol extract, exhibited strong cytotoxicity against A549 cells in a concentration-dependent manner, which was associated with apoptotic cell death as evidenced by morphological changes, Hoechst 33,258 nuclear staining, Annexin V and propidium iodide binding, and mitochondrial potentials assay 65 .New dammarane triterpene saponins called bacopaside É and bacopaside VII that were extracted from the n-BuOH fraction of Bacopa monniera demonstrated cytotoxicity against human tumour cell lines MDA-MB-231, SHG-44, HCT-8, A-549, and PC-3 M. Furthermore, at a concentration of 50 μmol/l, both compounds Table 3.The half-maximal inhibitory concentrations (IC 50 ) of the three pure compounds (1-3) isolated from Ziziphus spina-christi on four human cancer cell lines.GraphPad Prism 9 software was used to calculate the half maximal inhibitory concentration IC 50 (µg/mL) values.The concentration of each pure compound was transformed to log10.All experiments were performed at least two times.*p < 0.05, Student's t-test.www.nature.com/scientificreports/markedly reduced the adhesion, migration, and matrigel invasion of the breast cancer cell line MDA-MB-231 in vitro.These compounds at 50 μmol/kg demonstrated 90.52% and 84.13% inhibition in mice implanted with sarcoma S180 in in vivo studies 66 .
In human cancer cell lines, panaxadiol, a protopanaxadiol-type ginsenoside with a dammarane structure, also showed anticancer effects when combined with cyclophosphamide, 5-fluorouracil, and epicatechinin 73,74 .Strong anticancer activity has been observed for three derivatives of panaxadiol: 3β-acetoxy-panaxadiol, 3β-palmitic acid aceloxy-panaxadiol, and 3β-octadecanoic-panaxadiol 75 .High inhibitory effects against EBV-EA were demonstrated by dammarenolic acid, with IC 50 values of 226 mol ratio/32 pmol TPA.Ginsenosides are categorised as either dammarane-type or oleanane-type based on the differences in their aglycones.In the first group, it was demonstrated that 20(S)-ginsenoside Rg3 (Rg3) inhibited the proliferation of multiple tumour cells, including B16 melanoma cells, Lewis lung cancer cells, prostate cancer cells, and colon cancer cells 76 .
There are currently no known FDA-approved saponin-based anticancer medicines, despite the substantial research on saponins and their strong anticancer properties.This has a number of drawbacks, including toxicities and drug-like characteristics.In recent investigations, methods for ensuring enhanced efficacy and less toxicity in saponin were examined, including combination therapy and drug delivery systems 27 .

PPI network of the cancer related targets and KEGG-based enrichment analysis
The possible interactions between cancer-relevant proteins, first, all reported human-based proteins related to cancer disease, particularly for those tested in vitro in the present study (i.e., brain, breast, lung, colon cancers) were studied.Using each cancer type and "human"as keywords, we searched the Toxicogenomics (https:// ctdba se.org/) and GeneCards databases, as well as the previously published literature.As a result, more than 1043 proteins relevant to these cancers were retrieved.Only the common proteins (551 proteins) for each cancer type were selected for the subsequent analysis (Table S1).
We utilized the Cytoscape -software to create a protein-protein -interaction (PPI) network from these 551 selected proteins.
Figure 4A displays the network characteristics of the generated PPI, which included a high degree of connectivity (as seen by the 7602 edges between 306 nodes, a mean -node degree of 49.7, and a local clustering value of 0.557).The remaining selected proteins in 551 proteins list (Table S1) did not show any connections, so they were removed from the PPI network.
Therapeutic strategies for cancer may have a better chance of success if they focus on proteins having a high level of interaction, as these are typically the most significant and pertinent molecular targets (i.e., hub proteins or genes) in a given network 77 .Therefore, we focused on the top 8% (38 proteins) of the most intensively molecular targets that interact (i.e., hub proteins) sorted according to their degree value (Fig. 4B).
In addition, We organised the proteins in the existing network based on their role in the various signaling pathways associated with cancer pathogenesis.The KEGG database (https:// www.genome.jp/ kegg/ pathw ay.html) was used to guide this protein enrichment analysis.Proteins presented in the PPI network (Fig. 5) were categorized according to their involvement in the key cancer disease pathways into 9 groups: (i) Glioma signaling pathway; (ii) P3K-Akt signaling pathway; (iii) Adherens junction signaling pathway; (iv) Small cell lung cancer signaling pathway; (v) Autophagy signaling pathway; (vi) Ubiquitin mediated proteolysis pathway; (vii) Hippo signaling pathway; (viii) TGF signaling pathway; and (ix) TNF signaling pathway (Figs. 5, 6).Altogether, the current PPI network for human cancer offers a quick overview of the interacting proteins and associated signalling pathways, suggesting essential proteins that can be deemed critical to disease progression and thus suitable drug development targets.

Prediction of the target proteins for compounds 1-3
Multiple in silico-based experiments were subsequently performed on compounds 1-3 to putatively characterize its potential as an anticancer agent.The modeled structures of compounds 1-3 were run through the PharmMapper (http:// www.lilab-ecust.cn/ pharm mapper/) prediction platform to see if they might bind to any cancer-relevant protein target.PharmMapper is a unique pharmacophore-based virtual screening platform that can match the query structure into the 3D active sites-derived pharmacophore maps of most proteins hosted in the protein data bank (PDB; https:// www.rcsb.org/) 56 .
Unfortunately, none of the three compounds showed interactions with any cancer-related protein with an acceptable score (Fit score > 4).Even direct docking of the structures of these compounds to the cancer-relevant hub-proteins (Fig. 4) did not retrieve significant interaction with acceptable scores (< − 5 kcal/mol).These negative findings could be attributed to the large structure and the hydrophilic nature of these compounds that made them undruggable structures.The bioactivity of triterpenoidal and steroidal glycosides is usually linked to their aglycone part [78][79][80][81] because these glycosides can be easily hydrolyzed either spontaneously or in acidic conditions like this found in the stomach 82,83 .Hence, the aglycon parts of these glycosides are the key players in the bioactivity of such molecules e.g., glycyrrhetinic acid and diosgenin are the aglycones of both glycyrrhizin and dioscin, respectively.Both aglycones have been found to be the active anticancer molecules 79,84 .
Accordingly, we used the aglycone of the three compounds (1-3) for our in-silico analysis to link these major metabolites and their cancer-relevant target(s).So, the modeled structure of saponins aglycone was subjected to the PharmMapper-based pharmacophore screening to find out the potential cancer-relevant protein targets for this structure.Four proteins (i.e., BDR3, GLO1, CDK6, and HRAS) were predicted as potential targets for saponins aglycone only of which two proteins (i.e., CDK6, and HRAS) have been identified as hub proteins in the cancer PPI network described above (Fig. 4).Additionally, both proteins were involved in the PI3K-Akt signaling pathway (Fig. 5), indicating that saponins aglycone probably can kill cancer cells via this signaling pathway.
By generating another PPI network for the previously identified hub proteins (i.e., 38 highly interacting proteins in the network presented in Fig. 4) together with both GLO1 and BDR3, we were able to identify how the latter two proteins (i.e., BDR3 and GLO1) are essential in relation to these hub proteins.As shown in Fig. 7, BDR3 was found to interact with 25 hub proteins, while GLO1 did not show any interactions.This finding indicates that targeting BDR3 by saponins aglycone might also achieve a significant anticancer effect.
Saponin aglycones were then subjected to molecular docking and MD simulation experiments, which helped us to refine our preliminary pharmacophore-based virtual screening.Proteins were targets for saponins aglycone if they had docking scores with saponins aglycone below − 7 kcal/mol and absolute binding free energies (ΔG binding ) below − 7 kcal/mol.Accordingly, all the predicted proteins followed these conditions (Table S2), and hence, subjected to further refinement in 100 ns-long MD simulations.

Analysis saponins aglycone's modes of interaction
To investigate the binding modes of saponins aglycone inside the active site of each predicted protein, its modeled structure was re-docked into these active sites using Autodock Vina, and then the top-scoring binding pose for www.nature.com/scientificreports/ each protein was inspected and subjected to 100 ns-long MD simulation to check for the binding stability od saponins aglycone inside the active site of each protein.
First, BRD3: saponins aglycone and the co-crystallized ligand were found to be in a good alignment, and hence, were able to form stable binds within the enzyme's active site, with modest fluctuations and average RMSDs of 1.22 Å and 1.78 Å, as shown by 100 ns of MD simulations (Fig. 9A and B).Both saponins aglycone and the co-crystallized ligand shared the same hydrophobic interactions with TRP-332, PRO-333, VAL-338, LEU-343, LEU-345, and VAL-397, while differed in the H-bond formation, where saponins aglycone established single H-bond with the main chain of PRO-333, and the co-crystallized ligand formed 4 H-bonds with ASN-391.BDR3 (Bromodomain-containing protein 3) is a small globular protein that has been found to play a key role in chromatin remodeling and act as an epigenetic modifier 85,86 .Depletion or inhibition of BRD3 slows the growth of a number of cancer cell lines, including prostate and glioblastoma cell lines 87,88 .Recently, the withanolide physachenolide C (BDR3's co-crystalized ligand) has been reported to potently inhibit BDR3 89 .physachenolide C is a steroidal derivative, and hence, similar derivatives like saponins aglycone could be a good alternative inhibitor for this important cancer-relevant protein.
Second, GLO1: saponins aglycone and the co-crystallized ligand (i.e., 18-β-glycyrrhetinic acid; a triterpene derivative like saponins aglycone) were found to be in a good alignment.Accordingly, they were able to form similar interactions within the enzyme's active site (Fig. 9C).Both saponins aglycone and the co-crystallized ligand shared the same hydrophobic interactions with PHE-68, LYS-151, PHE-163, and TRP-171, while differed in the H-bond formation, where saponins aglycone established two H-bond with the LYS-151, and the co-crystallized ligand formed 3 H-bonds with ARG-38, ASN-104, and ARG-123 (Fig. 8C).Regarding their binding stability, both saponins aglycone and the co-crystallized inhibitor exhibited almost the same fluctuations and RMSD profiles (average RMSD = 2.34 Å and 2.84 Å, respectively) over the course of 100 ns of MD simulations (Fig. 9C).
GLO1 (glyoxalase I) catalyzes the conversion of hemimercaptal, formed from methylglyoxal and glutathione, to S-lactoylglutathione, and is involved in the regulation of TNF-induced transcriptional activity of NF-kappa-B 90 .GLO1 is overexpressed in cancer cells in response to the redox state, and hence, it is related to multi-drug resistance in chemotherapy, making GLOI modulators and/or inhibitors as potential anticancer therapeutics 91 .www.nature.com/scientificreports/Third, CDK6: saponins aglycone and the co-crystallized ligand (i.e., fisetin) were in a comparable alignment.Accordingly, they were able to form similar interactions within the enzyme's active site (Fig. 8C).However, the flavanol derivative fisetin (i.e., the co-crystalized inhibitor) established more H-bonds due to its multiple hydroxyl groups.Both saponins aglycone and the co-crystallized ligand shared the same hydrophobic interactions with ILE-19, VAL-27, VAL-77, PHE-98, and LEU-152, while differed in the H-bond formation, where saponins aglycone established two H-bond with the main chains of VAL-101 and GLN-103, and the polyhydroxylated co-crystallized inhibitor formed 8 H-bonds with LYS-43, GLU-61, GLU-99, ASP-163, and ASP-104 (Fig. 8C).About their binding stability upon simulation, both saponins aglycone and the co-crystallized inhibitor exhibited comparable fluctuations and RMSD profiles (average RMSD = 1.87 Å and 2.49 Å, respectively) over the course of 100 ns of MD simulations (Fig. 9C).
CDK6 (cyclin-dependent kinase 6) CDK6 is involved in a positive feedback loop that activates transcription factors through a reaction cascade.Importantly, these C-CDK complexes act as a kinase, phosphorylating and inactivating the protein of Rb and p-Rb related "pocket proteins" p107 and p130.While doing this, the CDK6 in conjunction with CDK4, act as a switch signal that first appears in G1, directing the cell towards S phase of the cell cycle 92,93 .CDK6 is overexpressed in cancer cells, and hence, it has been considered a key target for the development of efficient anticancer agents 94 .
Finally, HRAS: saponins aglycone and the co-crystallized ligand were not fully aligned with each other.(Fig. 8C).Hence, they shared only the same hydrophobic interactions with LEU-56 and TYR-71, while saponins aglycone did not form any H-bonds inside the active site of the enzyme (Fig. 8C).However, it was able to establish stable binding comparable to the co-crystallized inhibitor over the course of the MD simulation with an acceptable deviation (average RMSD = 2.66 Å and 1.99 Å, for saponins aglycone and the co-crystallized inhibitor, respectively) (Fig. 9D).
HRAS (Harvey rat sarcoma viral oncogene homolog) binds GDP/GTP and possesses intrinsic GTPase activity that is involved in the activation of Ras protein signal transduction 95 .This oncoprotein has been proven to be associated with a wide variety of tumors, and hence, targeting and inhibiting this protein is an excellent strategy to develop efficient anticancer medications.Until recently, there were not previously described HRAS inhibitors (similar to the remaining RAS oncoproteins).In 2020, several new potent HRAS inhibitors were reported targeting an allosteric binding site 96 , which in the present study, was also predicted to be a probable binding site for saponins aglycone.
From the previous network pharmacology-based and in silico-based findings, it can be concluded that saponins aglycone is a promising scaffold for the development of novel and efficient anticancer agents that can act via multiple targets and signaling pathways.Accordingly, such a comprehensive bioinformatics-based investigation should be the primary step in the development of new therapeutics, particularly from natural products.
On the other hand, it should be emphasized that such in-silico analysis presented herein has some limitations and assumptions that need to be considered when interpreting the results.First, the docking methods used to predict the binding poses and affinities of the compounds are based on empirical scoring functions that may not accurately capture the complex interactions and dynamics of the protein-ligand system.Second, the docking methods do not account for the effects of solvation, entropy, and conformational changes on the binding free energy.Third, the docking methods rely on the availability and quality of the protein structures, which may not reflect the native conformation or flexibility of the target.Therefore, the obtained in silico based-results should be regarded as tentative hypotheses that need to be validated by experimental methods.

Conclusion
The findings of this study demonstrated that three pure dammarane type saponins compounds which identified as christinin A, E, and F, and, derived from Z. spina-christi leaf extract, exhibited anti-cancer efficacy againt A549, U87, MDA-MB-231, and CT-26 cell lines.Notably, these compounds were strongly active at concentrations below 10 μg/mL, resulting in a 50% cytotoxic effect.Our results showed that all compounds exhibited low IC 50 values, indicating their potency at low concentrations across all four cancer cell lines.A PPI network analysis was used to identify 17 potential targets of saponins aglycone in cancer cells, and the binding modes of saponins aglycone to four of these targets (BRD3, GLO1, CDK6, and HRAS) were confirmed using docking and MD simulations.It was found that saponins aglycone could bind stably and similarly to the co-crystallized ligands of these proteins, and modulate their functions and affect their cancer-related pathways.As a result, anticancer effects such as cell cycle arrest, apoptosis, DNA damage, and inhibition of cell migration and invasion were induced by saponins aglycone.These observations alongside those obtained from the comprehensive network pharmacology-based and in silico-based analyses suggest that these compounds, particularly their aglycone, hold promise as potential candidates for cancer treatment.However, further investigations are necessary to elucidate the cellular mechanism of action for these compounds.Beside, in vivo studies to validate the cytotoxic effects of the isolated saponins, as well as an evaluation of their toxicity and potential side effects.
1 H and13 C NMR spectra of 2 were virtually superimposable to those of 1, with the addition of two signals attributable to carboxylic groups in the DEPT-Q spectrum (δ C 159.0, and 174.7), hence, the saponin was esterified with an oxaloyl moiety, reflecting the additional 71.98 mass units, and 2 degrees of unsaturation increase observed for 2 versus 1.The signals corresponding to the C-4ˋˋˋˋ of the trans-p-coumaroyl moiety (δ C 165.1) appeared downfield compared to those in 1 (δ C 159.0), while the signals corresponding to the CH-2ˋˋˋˋ, 6ˋˋˋˋ & CH-3ˋˋˋˋ, 5ˋˋˋˋ of the trans-p-coumaroyl moiety (δ H 7.56, 6.81δ C 130.4, 115.0) appeared upfield compared to those in 1 (δ H 7.59, 6.78, δ C 133.6, 115.9), and thereby indicated attachment of the oxaloyl moiety at C-4ˋˋˋˋ in trans-p-coumaroyl moiety in 2. Consequently, saponin 2 was identified as christinin F (Figs.1, 2).

Figure 3 .
Figure 3.The anticancer activities of the three pure compounds (1-3) isolated from Ziziphus spina-christi using the MTS assay on glioblastoma, breast cancer, lung, and colorectal carcinoma cells.Cancer cells lines A549 (a) U87 (b), MDA-MB-231 (c) and CT-26 (b) were tested.Cells (2 × 10 4 cells per well in 96 well plates) were treated with pure compounds as using the following concentrations (1, 5, 10 or 20 μg/mL) or a matching solvent control for 24 h at 37 °C with 5% CO 2 .Cells were stained with MTS reagents and cell survival was quantitated using the spectramax plate reader.The graphs represent the analysis of six replicates and show the percentage of survival of A549, U87, MDA-MB-231 and CT-26--treated cells relative to matching solvent-treated cells.The solvent used is methanol.All experiments were performed at least two times.*p < 0.05, Student's t-test. https://doi.org/10.1038/s41598-023-46841-2www.nature.com/scientificreports/

Figure 4 .
Figure 4. (A) Human cancer PPI network.This network comprises 306 nodes and 7602 edges with an average node degree of 49.7.The top-interacting nodes were colored by red (7.9%, 38 proteins of all interacting nodes, i.e., hub protein).(B) The top interacting-nodes (i.e., hub nodes arranged by their degree value).Green arrows represent the proteins predicted as probable targets for saponins aglycone investigated in the presented study.The thickness of the lines (i.e., edges) represents the degree of confidence (i.e., the strength of data support).

Figure 5 .Figure 6 .
Figure 5. Subnetworks derived from the complete cancer-relevant network in Fig. 4.These networks represent the first 4 key pathways involved in cancer pathogenesis according to the presented PPI networking in Fig. 4.The thickness of the lines (i.e., edges) represents the degree of confidence (i.e., the strength of data support).Nodes encircled with red line in the PI3K-Akt signaling pathway cluster are CDK6 and HRAS that were predicted to be potential targets for compound 3.

Figure 7 .Figure 8 .
Figure 7. (A) PPI network showing the interactions among the 38 hub proteins characterized in the complete cancer PPI network in Fig.4including both GLO1 and BDR3 that were identified as a probable target for saponins aglycone.Lemon-green-colored nodes represent proteins that show interactions with BDR3, while the red-colored ones are not.GLO1 (the blue-colored node) do not show any interactions.(B) The structures of saponins aglycone along with the structures of the four proteins predicted to be potential target for it.